Is Random Close Packing of Spheres Well Defined? 
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Despite its long history, there are many fundamental issues concerning random packings of spheres 
that remain elusive, including a precise definition of random close packing (RCP). We argue that the 
current picture of RCP cannot be made mathematically precise and support this conclusion via a 
molecular dynamics study of hard spheres using the Lubachevsky-Stillinger compression algorithm. 
We suggest that this impasse can be broken by introducing the new concept of a maximally random 
jammed state, which can be made precise. 
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Random packings of identical spheres have been stud- 
ied by biologists, materials scientists, engineers, chemists 
and physicists to understand the structure of living cells, 
liquids, granular media, glasses and amorphous solids, 
to mention but a few examples. The prevailing notion 
of random close packing (RCP) is that it is the maxi- 
mum density that a large, random collection of spheres 
can attain and that this density is a universal quantity. 
An anonymous author summarizes this traditional view 
as follows: "ball bearings and similar objects have been 
shaken, settled in oil, stuck with paint, kneaded inside 
rubber balloons - and all with no better result than (a 
packing fraction of) ... 0.636." |Q 

One aim of this paper is to reassess this commonly held 
view. First, we observe that there exists ample evidence 
in the literature (in the form of actual and computer 
experiments) to suggest strongly that the RCP state is 
ill-defined and, unfortunately, dependent on the protocol 
employed to produce the random packing as well as other 
system characteristics. In a classic experiment, Scott and 
Kilgour H obtained the RCP value 4> c « 0.637 by pour- 
ing ball bearings into a large container, vertically vibrat- 
ing the system for sufficiently long times to achieve max- 
imum densification, and extrapolating the measured vol- 
ume fractions to eliminate finite-size effects. Important 
dynamical parameters for this experiment include the 
pouring rate and both the amplitude and frequency of vi- 
bration. The key interactions are interparticle forces, in- 
cluding (ideally) repulsive hard-sphere interactions, fric- 
tion between the particles (which inhibits densification), 
and gravity. It is clear that the final volume fraction can 
depend sensitively on these system characteristics. In- 
deed, in a recent experimental study ||, it was shown 
that one can achieve denser (partially crystalline) pack- 
ings when the particles are poured at low rates into hor- 
izontally shaken containers. 

Computer algorithms can be used to generate and 
study idealized random packings, but the final states 
are clearly protocol-dependent. For example, a popular 
rate-dependent densification algorithm achieves 4> c be- 
tween 0.642 and 0.649, a Monte Carlo scheme || gives 



(j> c 1=3 0.68, and a "drop and roll" algorithm || yields 
(j> c w 0.60. It is noteworthy that in contrast to the last 
algorithm, the first two algorithms produce configura- 
tions in which either the majority or all of the particles 
are not in contact with one another. We are not aware 
of any algorithms that truly account for friction between 
the spheres. 

However, we suggest that the aforementioned incon- 
sistencies and deficiencies of RCP arise because it is an 
ill-defined state, explaining why, to this day, there is no 
theoretical determination of the RCP density. This is to 
be contrasted with the rigor that has been used very re- 
cently to prove that the densest possible packing fraction 
4> for identical spheres is n/ \/l8 « 0.7405, corresponding 
to the close-packed face-centered cubic (FCC) lattice or 
its stacking variants J?J . 

The term "close packed" implies that the spheres are in 
contact with one another with the highest possible coor- 
dination number on average. This is consistent with the 
view that RCP is the highest possible density that a ran- 
dom packing of close packed spheres can possess. How- 
ever, the terms "random" and "close packed" are at odds 
with one another. Increasing the degree of coordination, 
and thus, the bulk system density, comes at the expense 
of disorder. The precise proportion of each of these com- 
peting effects is arbitrary and therein lies the problem. 
In what follows, we supply quantitative evidence of the 
ill-defined nature of RCP via computer simulations, and 
we propose a new notion, that of a maximally random 
jammed state. 

A precise mathematical definition of the RCP state 
should apply to any statistically homogeneous and 
isotropic system of identical spheres (with specified in- 
teractions) in any space dimension d. Although we dis- 
card the term "close packed", we must retain the idea 
that the particles are in contact with one another, while 
maintaining the greatest generality. We say that a par- 
ticle (or a set of contacting particles) is jammed if it 
cannot be translated while fixing the positions of all of 
the other particles in the system. The system itself is 
jammed if each particle (and each set of contacting par- 
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tides) is jammed ||. This definition eliminates systems 
with "rattlers" (freely roaming caged particles) in the 
infinite- volume limit. We recognize that jammed struc- 
tures created via computer algorithms |j] or actual ex- 
periments will contain a very small concentration of such 
rattler particles, the precise concentration of which is 
protocol-dependent. Thus, in practice, one may wish 
to accommodate this this type of a jammed structure, 
although the ideal limit described above is the precise 
mathematical definition of a jammed state that we have 
in mind. Nevertheless, it should be emphasized that it is 
the overwhelming majority of spheres that compose the 
underlying "jammed" network that confers rigidity to the 
particle packing. 

Our definition of the maximally random jammed 
(MRJ) state is based on the minimization of an order 
parameter described below. The most challenging prob- 
lem is quantifying randomness or its antithesis: order. A 
many-particle system is completely characterized statis- 
tically by the JV-body probability density function P(r N ) 
associated with finding the system with configuration r^. 
Such complete information is never available and, in prac- 
tice, one must settle for reduced information. From this 
reduced information, one can extract a set of scalar or- 
der parameters ipx> ip2, ipm such that < tpi < 1, Vi, 
where corresponds to the absence of order (maximum 
disorder) and 1 corresponds to maximum order (absence 
of disorder). The set of order parameters that one selects 
is unavoidably subjective, given that there is no single and 
complete scalar measure of order in the system. 

However, within these necessary limitations, there is a 
systematic way to choose the best order parameters to be 
used in the objective function (the quantity to be mini- 
mized). The most general objective function consists of 
weighted combinations of order parameters. The set of 
all jammed states will define a certain region in the n- 
dimensional space of order parameters. In this region of 
jammed structures, the order parameters can divided up 
into two categories: those that share a common mini- 
mum and those that do not. The strategy is clear: retain 
those order parameters that share a common minimum 
and discard those that do not since they are conflicting 
measures of order. Moreover, since all of the parameters 
sharing a common minimum are essentially equivalent 
measures of order (there exists a jammed state in which 
all order parameters are minimized), choose from among 
these the one that is the most sensitive measure, which 
we will simply denote by tp. From a practical point of 
view, two order parameters that are positively correlated 
will share a common minimum. 

Consider all possible configurations of a d-dimensional 
system of identical spheres, with specified interactions, at 
a sphere volume fraction cp in the infinite- volume limit. 
For every (p, there will be a minimum and maximum value 
of the order parameter ip. By varying (p between zero 
and its maximum value (triangular lattice for d = 2 and 



FCC lattice for d = 3), the locus of such extrema de- 
fine upper and lower bounds within which all structures 
of identical spheres must lie. Figure 1 shows a schematic 
(not quantitative) plot of the order parameter versus vol- 
ume fraction. Note that at <j> = 0, the most disordered 
(tp = 0) configurations of sphere centers can be realized. 
As the packing fraction is increased, the hard-core in- 
teraction prevents access to the most random configu- 
rations of sphere centers (gray region). Thus the lower 
boundary of ip, representing the most disordered con- 
figurations, increases monotonically with (p. The upper 
boundary of ip corresponds to the most ordered struc- 
tures at each volume fraction, e.g. perfect open lattice 
structures (tp — 1). Of course, the details of the lower 
boundary will depend on the particular choice of ip. Nev- 
ertheless, the salient features of this plot are as follows: 

(1) all sphere structures must lie within the bounds and 

(2) the jammed structures are a special subset of the al- 
lowable structures |fTo|| . We define the MRJ state to be 
the one that minimizes ip among all statistically homoge- 
neous and isotropic jammed structures. 

To support the aforementioned arguments, we have 
carried out molecular dynamics simulations using sys- 
tems of 500 identical hard spheres with periodic bound- 
ary conditions. Starting from an equilibrium liquid con- 
figuration at a volume fraction of (p = 0.3, we compressed 
the system to a jammed state by the well-known method 
of Lubachevsky and Stillinger |£| which allows the diame- 
ter of the particles to grow linearly in time with a dimen- 
sionless rate T. Fig. 2a shows that the volume fraction of 
the final jammed states is inversely proportional to the 
compression rate T. A linear extrapolation of the data to 
the infinite compression rate limit yields (p ~ 0.64, which 
is close to the supposed RCP value reported by Scott and 
Kilgour. 

To quantify the order (disorder) in our jammed struc- 
tures, we have chosen to examine two important mea- 
sures of order: bond-orientational order and translational 
order [O. The first is obtainable in part from the pa- 
rameter Qq and the second is obtainable in part from the 
radial distribution function g(r) (e.g., from a scattering 
experiment). To each nearest- neighbor bond emanating 
from a sphere, one can associate the spherical harmonics 
Yi m {9, ip), using the bond angles as arguments. Then Qq 
is defined by |l2| 

\ m=— 6 / 

where Yg m denotes an average over all bonds. For a com- 
pletely disordered system in the infinite- volume limit, Qe 
equals zero, whereas Qq attains its maximum value for 
space-filling structures (Q^ cc ~ 0.575) in the perfect 
FCC crystal. Thus, Q§ provides a measure of FCC crys- 
tallite formation in the system. For convenience we nor- 
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malize the orientational order parameter Q = Qq/ Q§ cc 
by its value in the perfect FCC crystal. 

Scalar measures of translational order have not been 
well studied. For our purposes, we introduce a transla- 
tional order parameter T which measures the degree of 
spatial ordering, relative to the perfect FCC lattice at 
the same volume fraction. Specifically, 

where n, (for the system of interest) indicates the aver- 
age occupation number for the spherical shell of width a5 
located at a distance from a reference sphere that equals 
the ith nearest-neighbor separation for the open FCC 
lattice at that density, a is the first nearest-neighbor dis- 
tance for that FCC lattice, and Nc is the total number of 
shells (here we choose 8 = .196 and Nc — 7). Similarly, 
n ideai anc i n FCC are corresponding shell occupation 
numbers for an ideal gas (spatially uncorrelated spheres) 
and the open FCC crystal lattice. Observe that T = for 
an ideal gas (perfect randomness) and T = 1 for perfect 
FCC spatial ordering. 

The relationship between translational and bond- 
orientational ordering has heretofore not been character- 
ized. We have measured both T and Q for the jammed 
structures generated by the Lubachevsky-Stillinger algo- 
rithm and have plotted the results in the Q-T plane in 
Fig. 2b [^3|. This order plot reveals several key points. 
First, we observe that T and Q are positively-correlated 
and therefore are essentially equivalent measures of order 
for the jammed structures. Therefore, in seeking to deter- 
mine the MR J state using T and Q, one would search for 
jammed structures that minimize Q, the more sensitive of 
the two measures. Our preliminary results indicate that 
the MRJ packing fraction 4>mrj ~ 0.64 for 500 spheres 
using the Lubachevsky-Stillinger protocol. It should be 
noted, however, that a systematic study of other proto- 
cols may indeed find jammed states with a lower degree 
of order as measured by Q. Moreover, we notice that the 
degree of order increases monotonically with the jammed 
packing fraction [JllJ . These results demonstrate that the 
notion of RCP as the highest possible density that a ran- 
dom sphere packing can attain is ill-defined since one can 
achieve packings with arbitrarily small increases in vol- 
ume fraction at the expense of small increases in order. 

For purposes of comparison, we have included in the 
order plot of Fig. 2b results for the equilibrium hard- 
sphere system for densities along the liquid branch and 
densities along the crystal branch, ending at the maxi- 
mum close-packed FCC state jH. Interestingly, the equi- 
librium structures exhibit the same monotonicity prop- 
erties as the jammed structures, i.e., T increases with 
increasing Q and the degree of order increases with the 
packing fraction. Note that neither Q nor T are equal 
to unity along the equilibrium crystal branch because of 



thermal motion. 

To summarize, we have shown that the notion of RCP 
is not well-defined mathematically. To replace this idea, 
we have introduced a new concept: the maximally ran- 
dom jammed state, which can be defined precisely once 
an order parameter ip is chosen. This lays the math- 
ematical groundwork for studying randomness in dense 
packings of spheres and initiates the search for the MRJ 
state in a quantitative way not possible before. Never- 
theless, significant challenges remain. First, new and ef- 
ficient protocols (both experimental and computational) 
that generate jammed states must be developed. Second, 
since the characterization of randomness is in its infancy, 
the systematic investigation of better order parameters 
is crucial. 
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Fig. 1 Torquato et al. 




FIG. 1. A schematic plot of the order parameter if) versus volume fraction 4> for a system of identical spheres with prescribed 
interactions. All structures at a given packing fraction (j>, must lie between the upper and lower bounds (white region); gray 
region is inaccessible. The boundary containing the subset of jammed structures is shown. The jammed structures are shown to 
be one connected set; although, in general, they may exist as multiply disconnected. Point A represents the jammed structure 
with the lowest density and point B represents the densest ordered jammed structure (e.g., close-packed FCC or hexagonal 
lattice for d = 3, depending on the choice for V). The state which minimizes the order parameter ip is the maximally random 
jammed (MRJ) state. 
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FIG. 2. Molecular dynamics simulation results for the hard-sphere system, (a) The reciprocal compression rate 
versus the volume fraction (f> of the final jammed state of hard spheres using the molecular dynamics compression protocol of 
Lubachevsky and Stillinger |)|. The jammed state occurs when the diameters can no longer increase in time, the sphere collision 
rate diverges, and no further compression can be achieved after relaxing the configuration at the jammed volume fraction. Each 
point represents the average of 27 compressions, and the dashed line is a linear fit to the data, which yields (f> ~ 0.64 when 
r _1 = 0. (b) The Q-T plane for the hard-sphere system, where T and Q are translational and orientational order parameters, 
respectively. Shown are the average values for the jammed states of Fig. 2a (circles), as well as states along the equilibrium 
liquid (dotted) and crystal (dashed) branches. 
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